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Abstract 

We report on a ptychographical coherent diffractive imaging experiment in the water 
window with focused soft X-rays at 500 eV. An X-ray beam with high degree of 
coherence was selected for ptychography at the P04 beamline of the PETRA III syn¬ 
chrotron radiation source. We measured the beam coherence with the newly developed 
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non-redundant array method. A pinhole 2.6 /im in size selected the coherent part of 
the beam and was used for ptychographic measurements of a lithographically man¬ 
ufactured test sample and fossil diatom. The achieved resolution was 53 nm for the 
test sample and only limited by the size of the detector. The diatom was imaged at a 
resolution better than 90 nm. 

1. Introduction 

Imaging in the water window energy range between the absorption edges of carbon 
and oxygen at 284 eV and 532 eV yields a high chemical contrast in objects mainly 
composed of carbon and its aqueous components (Larabell & Nugent, 2010). The 
coherent X-ray diffractive imaging (CXDI) method (Miao et al.^ 1999) applied with 
third-generation synchrotron sources has proven to be a useful tool in structural anal¬ 
ysis on the nanoscale (Chapman &; Nugent, 2010; Mancuso et al.^ 2010a; Vartanyants 
& Yefanov, 2014). In CXDI no lenses are used and this in principle allows to over¬ 
come the resolution limitations of conventional lens microscopes. Reconstruction of an 
object from measured diffraction patterns requires solving the well-known phase prob¬ 
lem. Iterative phase retrieval techniques have been successfully employed to solve the 
phase problem (Fienup, 1982; Marchesini, 2007). As a limitation, the CXDI technique 
requires the sample to be isolated and fully illuminated by the coherent X-ray beam. 
In order to study extended objects with X-rays and to improve the uniqueness and 
convergence of the phase retrieval process, ptychographic coherent diffractive imag¬ 
ing (PCDI) was employed (Rodenburg et a/., 2007). PCDI involves scanning of the 
X-ray beam along the object up to a desired field of view. A certain overlap of the 
illuminated areas is crucial to succeed with the phase retrieval and the image recon¬ 
struction (Bunk et a/., 2008). Conventional CXDI and earlier iterative algorithms for 
PCDI required precise knowledge of the probe for example the X-ray beam intensity 
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profile incident on the object. A priori knowledge of the probe in ptychography is 
no longer necessary with algorithms that retrieve object and probe simultaneously 
(Thibault et a/., 2008). Moreover, ptychography has become an excellent tool to char¬ 
acterize optical elements such as pinholes (Giewekemeyer et a/., 20106), zone plates 
(Thibault et a/., 2008), focusing mirrors (Kewish et a/., 2010), compound refractive 
lenses (Schropp et a/., 2010), X-ray waveguides (Giewekemeyer et al.^ 2010a) and 
effects on the phase of the wave field in the focus (Dzhigaev et al.^ 2014). 

Third generation synchrotron sources provide intense and highly coherent X-rays 
which are widely used for coherent diffractive imaging on the nanoscale. The degree of 
coherence may not be constant along the cross section of the X-ray beam (Vartanyants 
& Robinson, 2003). We determined the coherence properties of the X-ray beam before 
performing our PCDI experiment and selected the most coherent part of the beam to 
avoid degradation of contrast in the diffraction patterns. Partial coherence may cause 
artifacts in the image reconstructions (Vartanyants & Robinson, 2001) and limits the 
spatial resolution. The most direct strategy to determine the spatial coherence of soft 
X-rays from synchrotron sources is to perform a set of Youngs double slit experiments 
with different separation between the slits (Chang et al.^ 2000; Paterson et a/., 2001). 
At the same time the approach of coherence characterization by non-redundant arrays 
(NRA) recently applied to X-rays (Skopintsev et al.^ 2014) offers a fast and reliable 
method to obtain the spatial coherence properties of synchrotron radiation from a 
single measurement. 

Fossil diatoms are unique monads with a light weight exoskeleton consisting of 
silicon dioxide (Si02). They exhibit very fine periodic three dimensional structures 
on the nano- and microscale at the same time which can hardly be manufactured by 
current nano-technology methods. A fossil diatom can thus be used as X-ray resolution 
test sample made by nature. Previously silica shells of fossil diatoms were successfully 
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studied with coherent imaging methods at synchrotron (Giewekemeyer et al.^ 20106) 
and free-electron laser sources (Mancuso et a/., 20106). The diatom investigated in 
this work belongs to the dominating species of nano-planktonic pennate Fragilariopsis 
cylindrus that is typically found in ice-edge zones in Antarctic waters (Kang &; Fryxell, 
1992). 

In this paper, we first present the characterization of the spatial coherence of the 
soft X-ray beamline with an NRA. This is followed by two PCDI measurements in the 
water window with optimized beam coherence. A lithographically manufactured test 
pattern of known structure and the fossil diatom are reconstructed as high resolution 
amplitude and phase contrast projection images. 

2. Experiment 

The experiment was performed at the soft X-ray beamline P04 (Viefhaus et al.^ 2013) 
at the PETRA III synchrotron radiation facility at DESY in Hamburg. The schematic 
layout of the beamline is shown in Figure 1(a). An APPLE-II type helical undulator of 
5 m length with 72 magnetic periods was tuned to deliver photons at an energy of 500 
eV which corresponds to a wavelength of A = 2.5 nm. The beam propagated to the 
dedicated X-ray vacuum scattering chamber [Holografische Roentgen-Streuapparatur, 
HORST, (Gorniak &; Rosenhahn, 2014)] through several optical elements, including a 
beam-defining slit (27 m downstream from the undulator), a horizontal plane mirror 
(35 m) and a monochromator unit consisting of a vertical plane mirror together with 
a plane varied-line-spacing (VLS) grating (46 m). The VLS grating focused the beam 
at the exit slit (71 m). A cylindrical mirror (79.1 m) collimated the beam in horizontal 
direction (not shown in Figure 1(a)). An elliptical mirror (78.5 m) focused the beam 
in vertical direction to the sample position (81 m). All mirrors were designed to accept 
a root mean square (rms) beam size of 6cr. 
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Knowledge about the coherence properties is an important prerequisite for CXDI 
experiments. For our ptychographic measurements we defined the size of the probe 
incident on the sample and selected the most coherent part of the X-rays at the same 
time with a pinhole of 2.6 /im diameter that was etched with a focused ion beam in a 
2 /im gold layer on a 100 nm thin Si 3 N 4 membrane. 


a 


HORST 


Undulator 



Beam-defining slit 


Plane mirror 

Exit slits Focussing mirror 



Plane VLS grating 


vacuum 

chamber 


b Coherence measurement 



Q Ptychography measurement 



Fig. 1. (a) Soft X-ray beamline layout. The X-ray scattering vacuum chamber HORST 

in (b) Coherence measurement setup and (c) ptychography setup. 

In both cases of the coherence and ptychography measurements the same sample 
holder of the HORST chamber was used (see Figure 1(b) and (c)). It consisted of a 
closed-loop piezo-electric stage to allow horizontal and vertical scanning of the sample 
relative to the probing X-ray beam with accuracy below 20 nm. The far-held diffrac¬ 
tion pattern intensities were measured by a CCD detector (DODX436-BN, Andor 
Technology Ltd., Belfast, UK). The square detector area of 27.6 x 27.6 mm^ consisted 
of 2048 X 2048 pixels with a pixel size of 13.5 x 13.5 /im^. 
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3. Coherence Characterization 

3.1. Theory 

A brief description of coherence theory is given in the following section to explain 
how the coherence properties of the beam can be retrieved from a single NRA diffrac¬ 
tion pattern (Skopintsev et al.^ 2014). In the theory of optical coherence, the statisti¬ 
cal properties of the radiation are described by the mutual coherence function (MCF) 
ri 2 (^) (Goodman, 2000; Mandel & Wolf, 1995) 

ri2{T) = {E*{ri,t)-E{r2,t + r)), ( 1 ) 

where and £’(r 2 ,t -h r) are the field values at positions and times ri,t and 

r 2 ,t -h r, and the brackets (• • •) indicate the average over time. The intensity li at 
position r^ is given by (|£’(r^, t)|^). The complex degree of coherence 7 i 2 (t) is defined 
as the normalized MCF 


712 (t) 


_ Fi2(r) _ 

V{|E(ri,t)R-(|E(r2,t)n' 


( 2 ) 


When the time delay r is much shorter than the coherence time Tc, the complex 
degree of coherence can be approximated by the complex coherence factor (CCF) 
712 = 712 ( 0 ) (Goodman, 2000). To characterize coherence by a single quantity the 
global degree of coherence is often introduced as (Vartanyants & Singer, 2010) 


. ^ /l7i2|V(ri)/(r2)dridr2 , . 

/I(ri)dri • //(r2)dr2 

In the frame of the Gaussian Schell-model (GSM), which in most cases provides 
sufficient physical description of the synchrotron radiation, the intensity profile and 
the CCF are both considered to be Gaussian functions (Mandel & Wolf, 1995). In this 
model the partially coherent beam is characterized by the standard deviation a of the 
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beam size and its transverse coherence length /coh- The coherence length is defined as 
the standard deviation of the modulus of the CCF | 7 i 2 |- In the frame of GSM, the 
global degree of coherence ( from equation (3) can be expressed as (Vartanyants & 
Singer, 2010) 


C — (^coh/^) 


4 + 


- 1/2 


(4) 


A Non-redundant array of apertures can be used to measure the CCF. It is shown 
(Mejia & Gonzalez, 2007; Skopintsev et a/., 2014) that for narrow bandwidth radiation 
the intensity / (q) of the far-held interference pattern as a function of the momentum 
transfer vector q observed in a diffraction experiment with N apertures is 


N 


I{q:)=Is{cO{Co + Y.C^^J 




(5) 


Here /s'(q) is the diffraction pattern of a single aperture. Individual aperture sepa¬ 
rations are denoted by = —^j,i and ai^j — —c^j^i are the relative phases. For the 
analysis of the diffraction pattern from equation (5) its Fourier transform is used 


/(Ax) = /s'(Ax) (g) 


N 


Co5(Ax) + ^C,j 





( 6 ) 


Here S(x) is the Dirac delta function. Is (Ax) is the Fourier transform of a single 

aperture diffraction intensity /s'(q) and the symbol (g) denotes the convolution. The 

N 

coefficient Cq is dehned as Co = Z] 4^, where C is the intensity incident on the i-th 

i=l 

aperture. The coefficients Cij are equal to the mutual intensity F^j(O) at r = 0 


Cij — F^j(O) 



(7) 


To each individual aperture separation dij corresponds a single peak in equation 
(6) with its height being equal to Cij. The intensities li^ Ij together with peak heights 
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Cjj are used to obtain the CCF values from equation (7) 



( 8 ) 


3.2. Results of the NR A eoherenee measurement 

The coherence properties of the P04 beamline were measured with a single NRA 
diffraction pattern for each set of beamline parameters and analyzed using the approach 
described in the previous section. The coherence was determined for 50 /im, 100 /im 
and 200 /im monochromator exit slit opening Dq^. From the monochromator resolving 
power of E/ AE = A/A A = 6 • 10^ the temporal coherence length /t = A • A/A A was 
estimated to be 15 /im. Geometric considerations showed that the maximum opti¬ 
cal path length difference in our experiment was / = 0.2 /im in the region used for 
the coherence analysis. This was much smaller than the temporal coherence length 
and confirmed our approximation of the complex degree of coherence by the CCF 
in the previous section. The spatial coherence of the X-ray beam was obtained by 
measuring the diffraction pattern produced by the NRA at the detector positioned 1 
m downstream (see Figure 1(b)). The NRA consisted of A = 6 identical rectangular 
apertures. Each aperture was 0.8 x 0.25 /im^ in size and manufactured according to a 
Golomb ruler (Lam & Sarwate, 1988). Our NRA was a Golomb ruler of order 6 with 
each separation between two individual apertures being unique (see inset in Figure 
2(b)). A background-corrected diffraction pattern of the NRA is shown in Figure 2(a). 
The Fourier transform I (Ax) of the measured diffraction pattern from the NRA is 
presented in Figure 2(b) and shows 31 well separated peaks. 


lUCr macros version 2.1.6: 2014/10/01 



9 



Fig. 2. (a) Diffraction pattern from the NRA diffraction measurement at 500 eV and 
Dqs = 100 /im. The white rectangle indicates the area used for the analysis, (b) The 
Fourier transform of the NRA diffraction pattern. Both images are displayed on a 
logarithmic scale. (Inset) Optical microscope image of the NRA and its aperture 
separations are shown in microns. 

To determine the CCF we analyzed the area (51 x 2001 pixels) shown as the white 
rectangle in Figure 2(a). In this region 51 line scans were used to determine the peak 
heights Cq and Cij. This area corresponds to the part of reciprocal space where the 
contribution of the high harmonics of X-ray radiation from the undulator is minimal 
(Skopintsev et al.^ 2014). 

The intensities li in the focus of the beam were determined by a beam profile scan 
with a pinhole of 1.5 /im diameter and are shown in Figure 3(a-c). For the 50 /im exit 
slit opening the intensity profile has a narrow peak. The side lobe at the right side of 
the profile was the result of diffraction from imperfections of the exit slit edges. At 
the 100 /im exit slit opening the side lobe has almost disappeared because a different 
section of the exit slit edge was illuminated. At the 200 /im exit slit opening a broad 
profile is observed without effects from the slit imperfections. The relative difference 
of the photon flux for each exit slit opening was measured by a photo diode. For the 
200 /im slit opening the highest flux was observed. In the case of 100 /im and 50 /im 
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the flux was reduced by a factor of 2.5 and 8, respectively. The flux was expected 
to be linearly dependent on the exit slit opening. We attributed the deviation to the 
imperfections and uncertainty of the exit slit positioning system at exit slit openings 
smaller than 100 /im. 


Des = 50 /im Des = 100 /im Des = 200 /im 



Fig. 3. Results of coherence measurement for three exit slit openings at 500 eV. 
Black dots indicate measured data and dashed lines represent Gaussian fits, (a)-(c) 
Intensity profiles measured with scans of a 1.5 /im pinhole, (d)-(f) Modulus of the 
CCF |7(Ax)|. The gray shaded area in (b) and (e) indicate the coherent part of the 
beam selected by the beam defining 2.6 /im pinhole. The rms values a of the beam 
size obtained from Gaussian fits, the coherence length /coh as well as the values of 
the global degree of coherence determined from equation (4) are also shown. For 
the Gaussian fits we used the data points up to 9 /im in (e) and up to 7 /im in (f). 


The CCF was obtained from intensities li and 51 sets of Cq, Cij. For each set 
of parameters Go, Cij the modulus of the CCF as a function of the NRA aperture 
separation Ax was retrieved using equation (8). In Figure 3(d)-(f) the averaged CCF 
I'fijl with error bars denoting the standard deviation is presented. 
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The spatial coherence length /coh was obtained by the Gaussian approximation 
exp (—Ax^/2/^qj^) of the averaged CCF |7i,j|. As expected, the coherence length and 
the global degree of coherence decreased nearly linear with the exit slit opening. The 
coherence length in horizontal direction was determined to be 12 /im for all exit slit 
openings (not shown here). 

The characterization of the coherence properties of the X-rays was performed in 
order to use an optimal exit slit opening for ptychography. Although a 50 /im exit 
slit opening provided the largest coherent fraction of the beam it also attenuated the 
beam strongly. Finally for ptychography measurements, the most coherent part of the 
beam provided by the 100 /im exit slit with the modulus of the CCF > 0.8 was 
selected by the 2.6 /im diameter pinhole (see shaded area in Figure 3(b) and (e)). 

4. Ptychography 

In far field ptychography the diffraction pattern amplitude Ai{qx^qy) in the detector 
plane is defined as the Fourier transform of the product of the probe P{x^y) and the 
projected object function Oi{x^y) 

A{(lx,qy)^^{P{x,y)-Oi{x,y)}. (9) 

Here x and y are the coordinates in the sample plane, qx and qy being momentum 
transfer values (i.e. the coordinates in reciprocal space) and T denoting the Fourier 
transform operator. The index i spans along the individual scan position of the probe 
in the object plane. From the ptychographical reconstruction in transmission geometry 
(see Figure 1(c)) one retrieves the complex object function 0(x, y') — exp(mA;A 2 ;(x, //)) 
which depends on the wave number k = 27r/A, the object thickness Az{x^y) and the 
complex refractive index n = 1 — 6 + i(3. Using the definition for n the object function 
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can also be written as 

Q(^x^y) = 


( 10 ) 


Here the first term with the absorption coefficient (3 represents the projected ampli¬ 
tude and the exponent of the second term with the refraction coefficient 5 is the 
relative phase shift A(p{x^y) = k6Az{x^y). Both terms contain the object specific 
response. First, we examined a strongly scattering tantalum (Ta) test object in the 
form of a Siemens star (ATN/XRESO-50HC, NTT-AT, Japan). The Siemens star was 
a lithographically manufactured sample with a smallest feature size of 50 nm. Second, 
we measured a fossil diatom skeleton dispersed on a silicon nitride (Si 3 N 4 ) membrane. 
With the 2.6 ym pinhole the incident beam was shaped to a highly coherent and cir¬ 
cular probe. Both samples were scanned on a rectangular scanning grid with a step 
size of 800 nm in both x and y direction perpendicular to the beam propagation axis. 
The detector was positioned in the far field at a distance of 58 cm downstream from 
the sample. 



Fig. 4. Diffraction images averaged over all positions i for the Siemens star (a, b) 
and for the Diatom (c, d). Measured (a, c) and reconstructed diffraction patterns 
(b, d) are separated by the white line. The diffraction patterns are displayed on a 
logarithmic scale. 
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We implemented the ePIE algorithm (Maiden & Rodenburg, 2009) to reconstruct 
both objects^. Initially the object function for all positions i was set to zero and 
the probe started with random values in amplitude and phase. The diffracted signal 
from the Siemens star was visible up to the edges of the detector. The reconstruction 
succeeded well over the whole diffraction plane (see Figure 4(a, b)). The diffraction 
signal from the diatom was comparably weak and reached the noise level already below 
the edges of the detector (see Figure 4(c, d)). As a consequence of the limited detector 
dynamic range, the signal above 7 /im“^ was dominated by noise. To reduce the 
noise contribution in the diatom diffraction patterns we subtracted the total signal 
corresponding to the mean noise level of each pixel plus four times the standard 
deviation from this mean that was determined from ten dark images. 

4 . 1 . Reconstruction of the probe function 

We exploited the strength of the ePIE algorithm and reconstructed the complex 
valued probe along with the object simultaneously. For both samples the probe was 
reconstructed and compared to a simulation in order to quantify the relative distance 
between the pinhole and the sample inside the HORST chamber. In the right column 
of Figure 5 we show the reconstructed probes from the experiment together with a 
calculated ideal probe. The white lines indicate the pinhole size. The ideal probe wave 
field originated as a plane wave from the 2.6 /im pinhole. In the middle column of 
Figure 5 we show the amplitudes retrieved by forward and backward Fresnel propa¬ 
gation of the complex valued wave fields. In the left column we show the beam cross 
sections in the pinhole plane that result from the propagations. The probe incident on 
the Siemens star was determined to be behind the pinhole at a distance of 0.08 mm. 

It contained a fine structure visible as concentric fringes that is typical for near-held 

^ Matlab (version 2013b) code was used for the entire analysis of the ptychography data. All calcula¬ 
tions were executed on a graphics processing unit (GPU). 
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diffraction from a pinhole with a Fresnel number 8.5 (see Figure 5(c)). We repeated 
the same procedure and found the diatom sample position at a distance of 0.81 mm 
behind the pinhole (see Figure 5(f)). The propagation of the diatom probe wave field 
to the position at z = 0 produced a well-shaped pinhole of 2.6 /im diameter and 
practically constant amplitude across the pinhole. 

We also compared our probe reconstructions with the simulation (see Figure 5(g-i)). 
The result of the simulation showed excellent agreement with the amplitude of the 
wave field distribution obtained from the Siemens star probe even at short distances 
behind the pinhole. The only difference was a slightly better contrast for the simulated 
propagation. A comparison with the wave field obtained from the diatom showed the 
same general behavior. However, it was lacking sharp features and high contrast which 
was both present in the Siemens star and simulated propagation. We attribute this 
effect to a lower signal produced by the diatom sample. 

The comparison of the wave fields at the position of the diatom {z — 0.81 mm) 
showed similar amplitude distribution with a beam size of about 1 /im full width 
half maximum (FWHM). Although this was considerably smaller than the pinhole 
(indicated as dashed-dotted black lines in Figure 5(b, e, h)) the probe was still large 
enough to provide a sufficient overlap between adjacent illumination positions for 
ptychography. 

One interesting result of this investigation is that one has to be careful if a pinhole 
is used in a ptychography experiment. The beam profile can be significantly different 
from the pinhole shape depending on the sample position behind that pinhole. In our 
experiment for example, at z = 0.35 mm the beam amplitude was practically equal to 
zero in the center of the beam, at the same time at the position of 0.21 mm behind the 
pinhole the beam had a sharp and narrow maximum with the size of 0.2 /im (FWHM) 
as well as strong shoulders on both sides (see Figure 5(b, h)). Good conditions for 
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our experiment were in principle the distances from 0.6 mm to 1.6 mm where the 
beam had a single pronounced peak and did not contain much amplitude in the wings 
outside of the geometrical pinhole region. A sample position (smaller than 0.1 mm) 
very close to the pinhole may also be beneficial because phase oscillations result in a 
structured probe that can improve the ptychographic reconstruction process in some 
cases (Maiden et a/., 2013; Quincy et a/., 2005). 




0.8 1 1.2 
z , mm 



Fig. 5. Ptychographic probe function reconstruction from (a-c) the Siemens star and 
(d-e) the diatom. For comparison (g-i) an ideal probe originating as plane wave 
from the 2.6 /im pinhole, (a, d, g) Amplitude distribution at the position of the 
pinhole; (b, e, h) distribution of the wave field amplitude downstream the pinhole; 
(c, f, i) amplitude distribution of the probe amplitude at the position of the sample 
(Siemens star (c), or diatom (f, i). Dashed white lines indicate the position of the 
sample relative to the pinhole at z = 0. White lines indicate the pinhole size. 


4-.2. Test pattern reconstruction 

The dataset from the Siemens star consisted of 132 diffraction patterns (12 hori¬ 
zontal X 11 vertical) that covered an area of 10.4 x 9.6 jiw? around the center of the 
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test pattern. The reconstructed Siemens star amplitude is shown in Figure 6(a). Black 
color denotes the strongly absorbing parts made of 270 nm thick substrate layers from 
SiC, Si 3 N 4 , Ru including the test pattern with 500 nm thick Ta with an expected 
transmission of Tmin = 4.6 • 10“^. In the white areas we expected a high transmission 
of Tmax = 0.48 that was defined by the substrate only. These areas had zero phase 
variation, as can be seen in Figure 6(b). In the low transmission parts a phase of 6.39 
rad was expected, however, undefined phases were observed. That could be explained 
by the phase value being close to 27r and small fluctuations causing complicated phase 
wrapping, both of which prevent a reliable and quantitative analysis (Giewekemeyer 
et al.^ 2011). 

Performing a series of reconstructions we observed in some cases a wrong number of 
lines (other than 36) in angular direction while the radius of the concentric patterns 
did not deviate from the expectation. This effect was also observed previously (Burdet 
et a/., 2014) and appears in case of not accurately known sample to detector distance. 
With our detailed probe function analysis we finally determined an accurate sample 
to detector distance of 578.5 ± 10 fim and obtained the excellent reconstruction shown 
in Figure 6. 

We used multiple angular line scans to plot the reconstructed amplitude contrast 
C — (^max “ ^min) / (^max + ^min) ^s a function of Spatial frequency (see Figure 
6(c)), where A^ax and Amin are the averaged maximum and minimum amplitudes of 
each line scan respectively. The error bars represent the standard deviation from the 
averaged contrast. The contrast decayed slowly from 0.87 at 1 period//im to 0.7 at 
the maximum spatial frequency of 9.4 periods//im, indicating a very good visibility 
for all spatial frequencies present in the reconstructed object. The maximum spatial 
frequency corresponds to a half period length of 53 nm and was equal to the pixel 
size of the reconstruction. Since the diffracted signal was cut at the detector edges in 
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reciprocal space (see Figure 4(a, b)) the obtained resolution was limited only by the 
detector size. 
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Fig. 6. Ptychographic reconstruction of the Siemens star test pattern, (a) Ampli¬ 
tude and (b) phase image, (c) Contrast C between high and low transmission as 
a function of spatial frequency from the angular scans denoted by red lines in the 
amplitude image. 
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4-5. Fossil diatom reconstruction 

In the case of the fossil diatom the dataset consisted of 119 diffraction patterns 
(17 horizontal x 7 vertical) covering an area of 14.4 x 6.4 As before in the 

case of the Siemens star reconstruction, we used our detailed probe function analysis 
to determine the accurate sample to detector distance and in this way avoided the 
problems of reconstruction of periodic structures. The reconstructed amplitude of the 
fossil diatom is shown in Figure 7(a). The reconstructed and unwrapped phase of the 
sample is shown in Figure 7(b). The color scheme displays the relative phase shift map 
A(f(x^y) between the substrate (shown in white) and the diatom. Assuming uniform 
density of Si02 (2.2 g/cm^) we determined the refraction coefficient 6 at 500 eV photon 
energy to be 5 = 1.39-10“^ (Henke et al.^ 1993). This allowed us to convert the relative 
phase shift A(p{x^y) between the substrate and the diatom to a map of projected 
material thickness by applying the relation Az{x^y) — A(f{x^y)/{k6). Diatoms of the 
Fragilariopsis cylindrus species have a cylindrically curved shape (see (Kang & Fryxell, 
1992)). To avoid an ambiguous thickness profile we show the projected thickness in 
Figure 7 up to a height of 500 nm. 

Ten equidistantly spaced ribs to be 1 ym period are visible in the reconstructed 
images. The length and width of each rib were estimated with 3 ym and 250 nm 
respectively. The fine structure that appeared in the form of a perforation and which 
is well pronounced in the amplitude image in Figure 7(a) had a period of 200 nm in 
vertical direction. All the discovered features of the fossil diatom that we investigated 
here are specific for the species Fragilariopsis cylindrus and are comparable with earlier 
studies (Kang & Fryxell, 1992; Bertilson et al.^ 2009). 

From line scans, which were extracted from the phase reconstruction, we determined 
the FWHM values of the fitted error functions, (see Figure 7(d)). Using the FWHM 
values, the resolution of our ptychographic diatom reconstruction was below 90 nm 
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and, consequently, 30 % better than in a previously published paper (Giewekemeyer 
et al.^ 2011). 
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Fig. 7. Ptychographic reconstruction of the amplitude (a) and the phase (b) of the 
fossil diatom, (c) Thickness map and (d) FWHM values of two error function fits 
along the black lines indicated in the phase reconstruction in (b). 
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5. Summary 

An experiment for high resolution ptychographical imaging of extended samples in 
the water window at 500 eV with highly coherent X-rays from the P04 beamlime 
at PETRA III has been presented. The spatial coherence was characterized with an 
efficient NRA method and high coherence of the X-ray beam was demonstrated. The 
global degree of coherence in the vertical direction varied from 10% to 73% depending 
on the setting of the exit slit of the monochromator. For ptychographic measurements 
an optimal exit slit size of 100 fim was used providing a coherence length of 4.1 /im 
and a degree of coherence of 35% at the sample position. With these settings the 
most coherent part of the beam was selected by the 2 /im pinhole that produced high 
contrast diffraction patterns in our ptychographic experiment. 

We obtained important knowledge about the fine features of the probe function by 
the propagation analysis that allowed us to determine the sample position inside of 
the HORST vacuum chamber with high precision. This turned out to be important 
in imaging periodic samples to avoid the artifacts caused by an inaccurately known 
sample to detector distance. The ptychographic reconstruction of the Siemens star test 
pattern was obtained up to 53 nm resolution and was limited only by the detector 
size. In the case of the weakly scattering fossil diatom we obtained the amplitude and 
phase of the transmission function quantitatively with a resolution better than 90 nm. 
Finally, we exploited the phase reconstruction and obtained a quantitative thickness 
map with details that allowed a comparison with the diatom species reported in the 
literature. 

Ptychography in the water window with its sensitivity for oxygen and carbon con¬ 
trast is especially promising to image hydrated biological cells. With the cryo-extension 
of HORST vacuum chamber its capability to image extended biological samples on the 
nanoscale will become feasible. We are also planning to extend our research to three- 
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dimensional (3D) ptychography with high spatial resolution (Dierolf et a/., 2010). This 
will allow us to determine the fine structure of diatoms on the nanoscale in 3D. 
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